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Abstract 

The relativistic equivalent of the Schrodinger equation for a two particle bound state 
having the total angular momentum S is written in the form of a Lorentz covariant set 
of equations {Pi + P2 + ^'^)^iPi,P2', P)xsiPi,P2)=P^^{pi,P2', P)xs{pi,P2) where the 
operators are the components of a 4- vector quasipotential. The solution of this set 
is a stationary function representing the distribution of spins and internal momenta in 
a reference frame where the momentum of the bound system is P^. 

The contribution of the operators to the bound state momentum is assumed 
to be the 4-momentum of a vacuum-like effective held entering the bound system as 
an independent component. It is shown that a state made of free quarks and of the 
effective held has dehnite mass and can be normalized like a single particle state. 

The generalization to the case of three or more particles is immediate. 
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It is a well known fact that the bound state problem has a fully satisfactory description 
in classical quantum mechanics, but does not have a similar one in relativistic field theory. 

In quantum mechanics the existence of a bound state is conditioned by the presence of a 
potential well which ensures the localizability of the wave function in the space of the relative 
coordinates. In this context the bound state wave function is stationary and normalizable. 

In relativistic field theories the bound state is assumed to be the result of the continuous 
exchange of quanta between the constituents as it is formally expressed by the Bethe-Salpeter 
method The formalism is fully relativistic, but the number of constituents in the inter- 
mediate states of the iterative solution is indefinite. 

It is almost obvious that the connection between the two schemes can be made only by 
assuming that the interacting potential in quantum mechanics is the effective, time averaged 
result of the continuous exchange of quanta. The difference between the two formalisms is 
clearly a difference of scale. One expects then for field theories to be more adequate for 
the treatment of high resolution processes occuring at high energy where the binding effects 
are negligible, while quantum mechanics to be suitable for the description of low resolution 
processes where the quantum fluctuations cannot be observed as such. It remains still to 
find a suitable approach for the bound state to be used in the treatment of low resolution 
relativistic processes, where the binding effects may be significant. 

The approach we propose here is suitable for this last case. It is based on a covariant 
generalization of the time independent Schrodinger equation for a bound state. This is 
achieved by replacing the eigenvalue equation for the Hamiltonian of the bound system by a 
set of four eigenvalue equations for a 4-vector operator representing the sum of the internal 
4-momenta and of a 4-vector interaction potential. The equations are written in momentum 
space and their solution are the distribution function of the internal spins and momenta. 
Lorentz covariance is manifestly satisfied and also is the mass shell constraint of the bound 
state momentum. 

In the following we refer specifically to QCD and to the meson case as a quark- ant iquark 
bound state. The generalization to three and more particles is briefly discussed. The colour 
indices shall be omitted for simplicity and the binding potential will be assumed to be white. 

The equations to be satisfied by the distribution function of the internal momenta 
'^s{PiiP2] P)Xs{PiiP2 of a meson with spin S and momentum P with = are: 



(pt + pi^ + n^)^ipi,p2;P)xsiPuP2) = P^-^{pi,P2]P)xs{pi,P2) /i = 0,l,2,3 (1) 

where 2 are the on mass shell momenta of the two quarks. The scalar function \l/(pi,p2; P) 
represents the distribution of the quark momenta and Xs{PiiP2) is an expression involving 
Dirac spinors, 7 matrices and quark momenta having suitable transformation properties. The 
operators fi^ behaving like the components of a 4-vector represent a relativistic generalization 
in momentum space of the interaction term and shall be called a generalized quasipotential. 
The generic form of the operators Vl^ can be written with the aid of the operators 
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pS=rf,2-^(Ol,2.)P^ (3) 

K2 = i^e-'''Pa7^^''H'''^ (4) 

K2 = ^^'''''pyi2'^i2 (5) 

representing the relativistic covariant generalizations of the position, momentum, spin and 
angular momentum in the rest frame of the meson. We have then: 

n>^ = — Uo + IpI', UiU + [p^', U2]+ + ^[v^ Vi]+ + z[v^, V2]+ (6) 

where V, U are scalar hermitian operators. We notice that Wo.1,2 are even at time reversal, 
while Vi,2 are odd. 

In the rest frame of the meson the equations (0) take the form: 



pf + mf + A/p! + mi+Wo U(Pi,P2;M)x5(Pi,P2) = M ^{p^p^; M)xs{pi,P2) (7) 



{Pl+P2 + [PlMlU + [P2,W2]+ + 2[Vi,Vi]+ +2[V2 V2] + ) ^{Pl,P2)Xs{Pl,P2) = 0. (8) 

If i7 = in this frame the solution of reads 

^{pi,P2;M) = 6^^\pi+p2) i>ip)xsip). (9) 

In addition, if Uq depends on pi — p2 only, in the nonrelativistic limit eq. (|^) reduces to the 
Fourier transform of the usual Schrodinger equation for a bound state 

-n 
P 

+ Uo)'>P{p)Xs{p) = {M -1711-1712) ip{p)xs{p), (10) 

where fj, = is the reduced mass and p = \{j>i — P2) is half of the relative quark 

momentum in the rest frame of the bound state. 

We notice also that the operators Ai and A2 are now replaced by the orbital angular 
momentum L of the relative motion 

Lk = - ie^jkP' (11) 

and the spin-spin and the spin-orbit couplings are ai^2 ■ ^1,2 and ai^2 ■ L where cTj are Pauli 
matrices. 

Then, just like in quantum mechanics one can write the orthogonality relation for the 
solutions of the eigenvalue equation (p!0|): 

d^P i^Aiip) i^li'ip) =Smm'- (12) 
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Equations (|T0|) and ( [T^ ) have been written for completeness and also for clarifying the 
relation between the set (|I|) and the Schrodinger equation for a bound state. 

Taking advantage of the relativistic covariance of ([^) and (0,^ one may solve them in 
the rest frame and write their solutions in any other frame simply by replacing any scalar 
product Pi ■ Pj by its Lorentz covariant expression —pj'^ ■ pj^ = -^{P ■ pi) [P ■ pj) — [pi ■ pj). 

We recall that equation (|^) with the solution (P) and Uq a relativistic scattering kernel 
0, and/or a suitable confining potential 0, HI, |@] represent the starting point for most 
relativistic quark models. In all these cases dynamics is restricted to the relative coordinate. 
Our approach is more general. The relative motion does not have a special significance 
and the quarks can be really treated as independent particles. This makes a generalization 
of the set (|l|) to the baryon case very easy by including the contribution of a third quark 
momentum and writing the suitable relativistic coupling of spins and angular momenta. 

The relation between the Bethe Salpeter formalism and the present approach also be- 
comes clear by looking at their solutions written as perturbative series with respect to the 
interaction term. In the BS formalism the solution writes with the aid of relativistic Green 
functions having poles at positive and negative energies. This makes the number of parti- 
cles in the intermediate states to be indefinite and causes a lot of trouble. In our case the 
perturbative series must be written for the relativistic invariant equation 

{P,pt + P.p'i + P^r]'^)^(pi,p2; P)xs{pi.P2) = MH{p,,P2; P)xs{puP2) (13) 

which has been derived from the set (|I|) and coincides with (|^) in the meson rest frame. 
Its Green function is Lorentz invariant and has only positive energy poles. This makes 
the perturbative series to resemble the series in nonrelativistic quantum mechanics and the 
number of particles to be the same in the external and intermediate states. 

We analyse now the way one can use the solutions of the set (|l|) to obtain a suitable 
representation of a meson in terms of independent quark operators. It is obvious that a meson 
is not an ensemble of free quarks with the momentum distribution '$s{pi,P2', P) because the 
sum of the free quark 4-momenta does not coincide with the bound state 4-momentum. In 
order to have a real dynamical representation of the bound state the contribution of the 
operators to the meson momentum must also be included. 

The solution we found to this problem is to write the contribution of the operators O'^ 
as the 4-momentum of an independent component of the meson existing besides the valence 
quarks p[. To this end we introduce the notation 

n>^^{pr,P2;P)xs{Pi,P2) = {P'^-Pt-P'i) I d'Q ^{puP2;Q) E u,M)^sVsM) 

S1,S2 

= [ d'Q ^ip^,p2; Q) E u,,{p,)VsVsM)- (14) 

where is a suitable tensor product of 7 matrices and quark momenta. The single meson 
state is written then as follows: 

|M,(P)) = / rfV 1^d'p2—d''Q 6^^\p, +P2 + Q-P) 
J ei 62 
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xip{puP2;Q) E ^sM)^sVsM) ^KQ) <(pi)&Up2)|0) (15) 



where a"^ and are quark and antiquark creation operators and u, v are free Dirac spinors. 
is a vacuum-like effective field carrying the momentum = —Pi — P2 defined in the 



relation (|T^). According to its definition the field $ represents the collective, time averaged 
effect of all the virtual excitations responsible for the binding and as such it must be compared 
with the "bag" in bag models ||^. We notice however that the bag has positive potential 
energy characterized by a constant volume density B, while the effective field $ carries the 
4-momentum whose values are a priori limited only by the quark and meson mass shell 



constraints, (see eq. (|1^).) The existence of a momentum Q besides the energy component 
Qo may be considered the consequence of the imperfect cancellation of the vector momenta 
in the virtual processes generating the binding and is in agreement with the opinion that 
the potential is the average effect of a continuous series of quantum fluctuations. 

The relativistic equivalent of eq.(|12]) is now the orthonormality relation for the single 
meson state ([T5|) . We calculate it in the general case by using the commutation relations for 
the free quark operators and assuming the following expression for the vacuum expectation 
value of the effective field: 



(0| $(Qi) $+(^2) |0) = fd'X e^ («-«^)'' ^" = ^ d^'\Q^ - Q2) (16) 

where Lq is the range of the binding forces and T is a time sensibly larger than the time 
basis involved in the definition of the effective field. 

We notice that S{Qo — Q'q) in (|16|) induces a cumbersome 6{E — E') in the expression of 
the norm. In order to avoid it by preserving the manifest Lorentz covariance of eq. ([T6|) we 
write 

1^ 5(Qo -Q'o) = ^[ dX, e^(^(^)-^(^'))Xo ^Er ^^^^.(m-mo^o ^ | ^mm' (17) 
{ M\P') \M{P) ) = 2E (27r)3 S(^\P - P') Smm' J (18) 



^0 ^0 
and get immediately: 



where 

J = / d'p, ^ d% ^ d'Q 5^'\p, +p, + Q- P)|(^(pi,P2; Q)? 

(Pl + mi -P2 + 1712 ^ \ ^ 

In the above relations we have implicitly assumed that M and M' are discrete eigenvalues 
of the equation (0) with \M-M'\To » 1 and hence the integral in (|17D vanishes if M 7^ M' . 
Notice the remarkable fact that J does not depend on the rather arbitrary time Tq. 



5 



In the same way as above it can be shown that the expression of the norm will contain 
the highly singular factor 5{a — a')6^{0) if fl vanishes in i.e. if = aP'^ where ck is a 
scalar. 

Relation (^) shows that the bound state function of the many particle system represent- 
ing the meson can be normalized like that of a single particle if the integral J' converges. 

We expect for the internal function (fi{pi,p2; P) Us^{pi)TsVs2{p2) which represents an 
equilibrium distribution of spins and momenta in the bound system to have a stationary 
counterpart in the coordinate space which must coincide in the nonrelativistic limit with 
the internal wave function of the meson. Proceeding like in the BS formalism we define the 
meson wave function in coordinate space 

(^(f , y, X, t; P) = {0\ ^(f , t)Ts'ilj{y, t)$(X, t) \MiP, E)) (20) 
where the single meson state is given by (|TBp, ip is the quark field 

t) = -i- Y: J d'p"^ (bt, v{p, s) e^^^'-^' + u{p, s) ^^^^t+wA^ (21) 



and 



<I>(X, t- Qo) = / d^ge-^0»*+'3^$(g, Qo). (22) 

After a straightforward calculation we get 

(p{x, y, X, t- P) = e-'^*+^^-^ ^(f , X- P) (23) 

where 



(27r)^ r d?pi d?p2 
L^J 2ei 262 



xTr {TsiPi + 7m)Ts'{-p2 + ma)!^)) ¥^(^1,^2; P - Pi - p^)eiP^-i--^)e'PMy~x) (34) 

describes the internal meson structure in coordinate space. 

At the first sight the expression (0) does not coincide with the usual internal wave 
function of a bound state because of its dependence on the coordinate X associated to the 
effective field which has no classical equivalent. However, in the nonrelativistic limit cj) writes 
as: 

(^(f , X- P) = 5^^\X - ^)0o(i^, m P) (25) 

if if depends on the relative momentum pi — p2 only. Here 0o is really the internal wave 
function in nonrelativistic quantum mechanics because it represents the Fourier transform 
of the nonrelativistic momentum distribution function. 

Concluding this paper we remark that the present approach provides a relativistic non- 
perturbative solution of the bound state problem. The main differences from the older 
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relativistic approaches are the covariant generahzation of the time independent Schrodinger 
equation accompanied by the introduction of a 4- vector quasipotential and the the definition 
of an effective field as an additional, independent constituent of the bound system besides 
the valence quarks. The 4-momentum of the effective field represents the contribution of the 
interaction term to the hadron momentum. The solutions of the dynamical equations (|l|) 
define equilibrium distributions of quark spins and momenta inside hadrons and hence may 
be directly used to calculate electroweak formfactors in the soft limit where binding effects 
are important. 
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